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SUMMARY 


This  report  covera  the  firat  year  of  e  two  year  project  dealing  with  the 
development  of  an  OBOGS  computer  model  to  aid  in  design  and  apecif ication  of 
these  oxygen  generation  ayetema  for  the  Army' a  in-flight,  medivac,  and  field 
hospital  use.  A  model  which  predicts  the  performance  of  a  two  bed  system  for 
°2  -  N2  separation  has  been  delivered  and  installed  on  the  U.S.  Army 
Aeromedical  Research  Laboratory  (USAARL)  VAX  computer.  A  model  to  predict  the 
propagation  of  contaminants  in  these  beds  is  currently  under  development. 


STATEMENT  OF  PROBLEM 


The  ultlnate  coal  of  this  reaearch  ia  to  Inaura  proper  design  of  nolecular 
sieve  oxygen  generation  systans  for  the  U.S.  Amy's  in-flight,  nedlvac,  and 
field  hospital  use.  Specifically  the  reaearch  lnvolvee  further  developnent  of 
an  0B06S  nodal  to  include  the  affects  of  conteaiinanta  in  the  feed  air.  This 
OBOGS  nodal  can  be  used  to  optinize  and  design  OBOGS  systens  with  respect  to 
systen  paraneters  such  as  cycle  tine  and  bed  and  valve  dinenslons. 


MCKSBPIfflB 


The  elimination  of  the  need  for  liquid  oxygen  storage  for  in-flight  and 
portable  ground  use  by  means  of  oxygen  generation  equipment  ie  an  especially 
attractive  goal  to  the  military.  Such  an  accomplishment  would  not  only 
increase  safety,  but  would  relieve  the  logistics  requirements  of  LOX  supply 
which  is  costly  and  can  restrict  mission  durations.  Although  other  means  are 
available,  the  molecular  sieve  method  of  generating  enriched  oxygen  breathing 
gas  is  the  moet  feasible  for  in-flight  usage.  The  molecular  sieve  system 
weighs  little,  and  has  the  potential  of  high  reliability. 

A  molecular  sieve  system  utilizes  the  procese  of  physical  adsorption  to 
separate  air  into  its  components.  Although  a  detailed  discussion  of  adsorptive 
phenomena  is  beyond  the  scope  of  this  report,  a  brief  overview  of  the  procese 
is  in  order. 

Adsorption  is  the  process  by  which  particular  types  of  molecules  contained 
in  a  liquid  or  gas  are  selectively  collected  onto  a  solid  surface.  The  term 
"adsorbent”  refers  to  the  collecting  material  while  the  "adsorbate”  is  that 
material  which  is  collected.  In  physical  adsorption  (as  opposed  to  chemical 
adsorption),  Vender  Waal's  forces  create  the  bonds  between  adsorbent  and 
adsorbate;  thus,  the  interaction  is  weak.  The  solid  adeorbent  is  typically  a 
highly  porous,  granular  material  with  interior  pores  of  diameters  on  the  order 
of  . 4nm,  and  with  the  total  pore  volume  approaching  SOS  of  that  of  the  whole 
particle.  Clearly,  abundant  internal  storage  space  is  available  in  the  sieve 
material . 

Adsorptive  separation  is  the  process  by  which  the  molecular  sieve  system 
would  separate  air  into  its  components  and  thereby  produce  an  enriched  oxygen 
product  gas. 

To  maintain  continuous  operation  from  an  adsorption  plant,  molecular  sieve 
beds  must  be  periodically  "regenerated”  by  a  period  of  desorption  in  which  the 
adsorbate  is  removed  from  the  sieve  material.  There  are  several  ways  in  which 


desorption  may  be  accomplished.  "Thermal  Swine"  involves  heating  the  bed  to  a 
temperature  at  which  the  adsorptive  capacity  of  the  bed  material  Is  so  low  that 
the  adsorbate  must  leave  the  molecular  sieve  internal  surfaces.  "Pressure 
swing”  creates  the  same  effect  by  reducing  the  adsorptive  capacity  at 
essentially  constant  temperature.  "Inert  purge  stripping"  removes  the 
adsorbate  at  constant  temperature  and  pressure  by  the  passage  of  fluid 
containing  no  adsorbate  molecules,  in  which  the  adsorbate  is  soluble. 
"Displacement  desorption"  occurs  when  a  fluid  containing  a  high  concentration 
of  a  more  strongly  adsorbed  molecule  is  passed  through  the  bed.  The  pressure 
swing  adsorptive  cycle,  PSA,  is  of  primary  concern  in  this  work. 

Two  important  technological  developments  of  recent  decades  have  created  the 
possibility  of  in-flight  and  portable  air  separation  systems.  Firstly,  the 
domestic  development  and  production  of  synthetic  zeolites,  a  molecular  sieve 
material,  affords  the  availability  of  reproducible  adsorbents.  Secondly,  the 
inception  and  development  of  the  pressure  swing  adsorption  cycle  (PSA)  has  been 
found  to  yield  a  satisfactory  production  system. 

The  PSA  process  usually  employs  two  or  more  molecular  sieve  beds  and  cycles 
them  in  a  manner  such  that  a  continuous  flow  of  enriched  oxygen  product  gas  is 
attained.  Sufficient  gas  is  produced  to  meet  demand  and  to  regenerate 
exhausted  beds  via  pressure  swing  method.  By  use  of  an  appropriate  zeolite, 
oxygen  or  nitrogen  (for  fuel  tank  inerting)  can  be  concentrated  to  almost  any 
value  up  to  95*  for  02and  99.8*  for  M2*  The  product  gas  can  ba  delivered 
at  normal  room  temperature  and  at  some  nominal  pressure  as  desired.  The 
process  cost  is  basically  the  cost  of  air  compression  for  the  feed  gas. 

Traditionally,  PSA  systems  have  found  applications  in  industry  and  have 
been  very  large  systems  which  were  relatively  Insensitive  to  the  various  rate 
Phenomena  of  the  adsorption  process.  In  the  development  of  a  scaled-down 
system  for  portable  applications,  knowledge  of  transient  behavior  becomes 
exceedingly  Important  due  to  the  vastly  increased  sensitivity  of  system  output 
to  rate  phenomena  and  the  requirements  of  minimal  size  and  weight. 
Additionally,  the  effects  of  potentially  harmful  contaminants  and  chemical 
agents  which  might  be  encountered  during  usage  must  be  better  understood. 


CONTAMINANT  MODEL 

As  a  first  step  in  the  Modeling  of  an  OBOGS  system,  an  experimentally 
verified  model  of  a  single  bed  system  containing  a  zeolite  molecular  sieve  with 
contaminants  in  the  inlet  air  is  being  developed.  A  model  that  has  been 
previously  developed  for  two  component  02  -  N2  gas  mixtures  is  described  in 
appendix  A.  Preliminary  work  with  02  -  N2  -  CO  and  02  -N2  -C02 
mixtures  indicates  that  a  similar  model  with  a  third  component  might  be 
effective  for  contaminant  studies.  USAARL  is  presently  conducting  single  bed 
experiments  with  contaminants  such  as  C02,  CO,  C2H60  (ethanol),  NH3 
(ammonia),  C5H12  (pentane),  CH3Br  (methyl  bromide),  N02,  N2O4,  03 
(ozone),  CH4  (methane),  CC12F2  (freon  12),  acrolein,  oil  breakdown 
products,  and  other  battlefield  products.  These  experiments  are  essentially 
column  breakthrough  experiments  which  will  be  used  to  validate  the  computer 
model . 

Once  the  dynamics  of  a  single  bed  are  correctly  modeled,  an  overall  system 
model  for  a  two  bed  system  with  valves  and  mixing  plenum  will  be  developed. 
This  modeling  has  been  done  for  02  -  N2  mixtures.  A  description  is 
contained  in  appendix  B.  A  similar  model  is  being  developed  for  contaminant 


A  computer  program  of  the  OBOGS  model  described  In  appendix  B  has  been 
delivered  and  Installed  on  the  VAX*  computer  at  USAARL.  This  program  is  menu 
driven  for  ease  of  operator  use  and  has  REGIS*  graphical  output  for  use  with 
VT125*  or  VT240*  graphical  terminals.  The  operator  has  a  choice  of  output  of 
on-line  plotting  of  oxygen  mole  fraction  or  Inlet/outlet  mass  flowrate  versus 
time  or  off-line  plotting  of  a  dynamic  simulation  of  OBOGS  for  one  cycle.  The 
operator  can  also  select  and  change  a  variety  of  geometric  and  material 
parameters  in  order  to  aid  in  design  studies.  Sample  inputs  and  outputs  of  the 
program  are  given  in  appendix  C. 

A  preliminary  3  component  OBOGS  has  been  developed  for  contaminant 
studies.  Before  any  further  development  on  this  program  can  continue  the 
results  of  the  experimental  single  bed  data  generated  by  USAARL  is  necessary. 
As  soon  as  this  data  is  available,  it  will  be  incorporated  into  a  single  bed 
model.  Then  a  two  bed  OBOGS  system  model  will  be  developed  and  experimentally 
verified. 
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APPENDIX  A 


Single  Bed  Model  Technical  Description 
This  article  is  reprinted  from  the  TRANSACTIONS  OF  THE  ASME,  JOURNAL  OF 
DYNAMIC  SYSTEMS,  MEASUREMENT,  AND  CONTROL.  The  authors  are  J.J.  Beaman, 
A . J .  Healey,  and  J.  Werlin,  Dec.  1983 


A  Dynamic  Model  of  a  Molecular 
Sieve  Bed  With  Nonlinear  and 
Coupled  Isotherms1 

A  mathematical  model  of  the  0}-N2  separation  process  in  a  molecular  sieve  bed  was 
developed  and  experimentally  verified.  Oxygen  and  nitrogen  breakthrough  ex¬ 
periments  were  performed  on  a  single  bed  packed  with  molecular  sieve  5 A  .  In  these 
experiments,  the  bed  had  initially  adsorbed  either  pure  oxygen  or  nitrogen  and  then 
a  constant  inlet  flow  of  nitrogen  or  oxygen  was  switched  on  the  bed.  The  outlet 
oxygen  mole  fraction  and  mass  flowrate  were  then  recorded.  This  was  done  at 
several  pressures  and  flowrates.  The  model  predictions  agreed  with  the  ex¬ 
perimental  data  when  a  linear  oxygen  isotherm  and  a  nonlinear  Langmuir  isotherm 
for  nitrogen  were  used  and  when  these  were  coupled.  The  numerical  scheme 
presented  in  this  study  to  solve  the  system  equations  makes  efficient  use  of  com¬ 
putation  time. 


I  Introduction 

The  elimination  of  the  need  for  liquid  oxygen  for  crew 
breathing  requirements  in  aircraft  is  an  especially  attractive 
goal  A  means  of  accomplishing  this  goal  is  now  available 
through  the  use  of  a  molecular  sieve.  OBOGS  (On-Board 
Oxvgen  Generation  System).  Molecular  sieves  are  in 
widespread  use  industrially  |l),  but  the  special  problems 
associated  with  an  OBOGS  have  led  to  the  need  for  a  reliable 
computer-based  simulation  tool  that  would  model  the 
re'atively  rapid  transient  response  of  such  systems. 

One  promising  system  uses  a  synthetic  zeolite  to  adsorb 
nitrogen  from  engine  bleed  air,  thereby  producing  oxygen 
enriched  breathing  gas.  This  system  uses  a  “pressure  swing” 
adsorption  cycle  [2|  in  order  to  achieve  gas  separation.  In  this 
cycle,  nitrogen  is  adsorbed  from  engine  bleed  air  at  high 
pressure  (producing  oxygen  enrichment)  and  then  desorbed 
and  exhausted  to  the  atmosphere  at  low  pressure  in  order  to 
replenish  the  sieve  material.  While  large  scale  industrial 
systems  are  in  commercial  use.  size  and  weight  optimized 
units  for  airborne  applications  are  still  in  the  development 
stage 

In  order  to  better  understand  the  operation  and  scaling  laws 
of  the  OBOGS.  a  computer  model  of  this  system  has  been 
developed.  As  a  first  step  in  this  development,  a  model  of  a 
single  molecular  sieve  bed  with  identified  parameters  is 
needed.  This  paper  addresses  the  formulation  and  verification 
of  such  a  model  for  predicting  the  transient  response  of  gas 
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flow  rate  and  gas  concentration  at  the  bed  outlet  to  sudden 
changes  in  gas  composition  at  the  bed  inlet. 

Previous  work  in  this  area  has  treated  the  diffusion  process 
to  be  instantaneous  [3,  4)  or  has  considered  one  component  at 
the  trace  level  [5).  This  work  considers  a  complete  “washout" 
of  one  gas  type  using  another.  In  other  words,  mole  fractions 
at  the  inlet  undergo  step  changes  from  0  to  I  and  vice  versa, 
while  the  bed  back-pressure  and  inlet  flows  are  essentially 
constant. 


II  Theoretical  Model 

Figure  1  shows  a  schematic  of  the  experimental  single  bed 
which  was  used  to  validate  the  molecular  sieve  bed  model, 
where  component  A  is  oxygen.  B  is  nitrogen,  jr  is  oxygen  mole 
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steve  fed  i-  25  4  cm  *  |4.0<m  x  14.5  cm.  loaded  with  3  8  kg 
oi  l  ,r.je  5A  sorbeni  screened  through  a  20  to  40  meshes  in 
This  corresponds  to  zeolite  particles  340  to  833  microns  in 
diameter  The  void  fraction  of  the  bed  is  approximately  <  = 
0.36  \I1  experiments  were  run  at  room  temperature  of  25  "C 

On  the  supply  side  of  the  bed  were  cylinders  of  pressurized 
O;  and  N;  Inflow  was  regulated  by  standard  reducing  valves 
and  two  variable  orifice  needle  valves,  each  of  which  con¬ 
trolled  the  flowrate  out  of  one  cylinder  only.  A  three-way  tap 
valve  was  used  to  produce  the  step  change  in  gas  composition. 
Pressure  transducers  were  located  at  the  bed  inlet  and  outlet. 
These  transducers  were  located  at  the  bed  inlet  and  outlet. 
These  transducers  served  as  a  check  on  the  isobahc  assump¬ 
tion.  The  amount  of  pressure  drop  across_the  bed  was  on  the 
order  of  3  percent  of  inlet  pressure.  A  mass  flowmeter  was 
located  lust  beyond  the  bed  outlet  to  measure  outflow.  Outlet 
0;  mole  fraction  was  measured  with  a  respiratory  mass 
spectrometer.  Experimental  data  were  taken  for  various 
flowrates  at  four  different  outlet  pressures:  W  atm,  1  atm,  2 
atm.  and  3  atm.  A  hypobaric  chamber  was  employed  for  the 

atm  test,  while  a  hyperbaric  chamber  served  as  the  plenum 
chamber  for  the  2  and  3  atm  tests.  For  more  information  on 
the  test s  see  Werlin  [17). 

Figure  3  shows  the  general  nature  of  the  oxygen  and 
nitrogen  ’  breakthrough”  experiments  where  the  symbols 
refer  io  experimental  data.  From  0  to  16  seconds  an  oxygen 
breakthrough  is  shown.  As  seen  in  this  figure,  the  outlet  mass 
flowrate  rises  to  approximately  twice  the  inlet  flowrate  and 
slowly  trails  off  to  its  equilibrium  value.  The  outlet  mole 
fraction  concentration  starts  at  pure  nitrogen  and  rises  slowly 
to  pure  oxygen  At  16  seconds,  the  inlet  flow  is  suddenly 
switched  to  pure  nitrogen.  As  this  point  the  outlet  flow  drops 
almost  immediately,  is  constant  for  a  while,  and,  finally  it 
rises  back  to  its  equilibrium  value  as  the  concentration 
"Iront”  propagates  through  the  bed.  The  concentration  time 
profile  is  much  sharper  than  the  oxygen  breakthrough  case, 
and  at  approximately  24  seconds  the  oxygen  has  been  com¬ 
pletely  washed  out. 

This  general  behavior  can  be  explained  by  the  fact  that  5A 
sorbent  has  a  larger  affinity  for  nitrogen  than  oxygen,  so, 
when  oxygen  is  washed  out,  more  nitrogen  is  adsorbed  than 
oxygen  is  given  up.  Therefore,  the  flowrate  will  be  reduced 
ahead  of  the  concentration  front.  The  concentration  front 
remains  sharp  due  to  the  fact  that  the  velocity  ahead  of  the 
front  is  lower  than  that  behind  it. 

For  oxygen  breakthrough,  the  effect  is  reversed.  In  this 
case,  more  nitrogen  is  desorbed  than  oxygen  is  adsorbed. 
There  is  a  higher  flowrate  ahead  of  the  front,  and  oxygen  and 
nitrogen  will  tend  to  be  convected  away  from  this  region, 
leading  io  a  broader  concentration  front. 

IV  Identification  of  Model  Parameters  and  Com¬ 
parison  With  Experimental  Data 

Effect  of  Model  Parameters.  A  quantitative  assessment  of 
the  efiect  of  the  model  parameters  on  the  characteristic 
flowrate  and  mole  fraction  time  plots  can  be  made  by  writing 
(5)  through  (8)  in  nondimensional  form  as 
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u  *  u/Uo,i*o  *  inlet  convective  velocity,  f  =  3Z/u0, 
ij  =  n/n0,  n,  »  initial  adsorbed  phase  concentration, 

1/t  *  nA^n i.  r  *  tKAd,  x  —  CA/C,  g  =  y-C,  and 

8  a  “  IM  fC. 

The  single  component  isotherms  of  (18)  and  (19)  can  be  ex 
pressed  as 

*  n°A  (28 a) 

8°  -  (I  -|i)t|J/(1  -  anj®)  (286) 

where  n  »  1  -  KflKA  and  a  *  n0/b  The  ideal  adsorption 
method  is  then  used  to  couple  these  isotherms.  For  a  nitrogen 
breakthrough,  the  initial  and  boundary  conditions  become 

w(0,t)  *  1 
xfO.r)  «  0 

■Hf.O)  =  I 

iiU.0)  =  1 

Equations  (24)  through  (28)  with  these  conditions  can  then 
be  solved  numerically,  as  discussed  previously,  for  values  of 
the  two  parameters  m  and  a  These  simulations  indicate  that 
the  nondimensional  distance  down  the  bed  is  related  to 
nondimensional  travel  time  by 

f«(l -x*oMr,0  -0  70)  (29) 

where  is  the  nondimensional  time  for  the  outlet  oxygen 
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Table  1  Nitrogen  breakthrough  experimental  results 


p 

Bed  pressure 
MPx 

«o 

Inlet  flowrate 
is 

w* 

Mass  flow  droop 

<50 

lime  to  50^1 
s 

0.0507 

3.21 

0.41 

9.40 

0  0507 

4.28 

0.43 

7.35 

o  too 

4.28 

0.50 

10.7 

0  100 

8.40 

0.52 

5.80 

0  204 

6  53 

0.63 

15.0 

0  203 

8  46 

0.65 

12.2 

0.303 

6.42 

0.70 

19.5 

0  303 

8.56 

0.72 

14.9 

mole  fraction  to  drop  to  0.5  or  50  percent.  This  linear  relation 
arises  due  to  the  fact  that  the  outlet  velocity  drops  to  a  value 

Uj  a  1  -  n  +  a  for  0  <  H  -  or  <  1 ,  (30) 

where  ud  is  the  droop  velocity  as  determined  by  simulation 
and  seen  in  Fig.  3  for  I  >  16  seconds.  This  velocity  remains 
relatively  constant  for  a  time  determined  by  the  length  of  the 
bed.  (The  0.70  in  the  (30)  above  comes  from  an  intrinsic 
exponential  time  constant  in  the  process,  i.e.,  e~  70  =  50 
percent.)  For  a  given  value  of  rM,  the  length  at  which  the  mole 
fraction  “front"  is  50  percent  would  be  approximately:  length 
=  droop  velocity  x  time  to  50  percent  which  is  the  in¬ 
terpretation  of  equation  (29). 

Identification  of  Model  Parameters.  The  values  of  KA , 
A*,  b,  and  d  can  be  effectively  found  by  using  breakthrough 
data  and  the  results  of  the  preceding  section.  The  steps  to 
obtain  these  parameters  from  a  nitrogen  breakthrough  are: 

( 1 )  Obtain  values  of  p  -  a  from  flowrate  droop  data,  i.e., 
*1-0=1  -  ud  for  at  least  two  different  bed 
pressures.  This  allows  calculation  of  *i  and  a. 

(2)  Use  (29)  to  obtain  a  value  of  KA  from  the  time  for  50 
percent  breakthrough.  The  value  of  d  is  superfluous  to 
the  calculation  since  0  »  d(l  -  »)/«  and  r^  >  >  0.70 
for  the  system  under  study.  Using  these  two  relations, 
(29)  yields 

A„  -L(  1  -<)/m0(1  ~  n  +  ahtio 

(3)  UseAs  =  (I  -  h)Ka  to  calculate  A*. 

(4)  From  dimensional  analysis,  one  sees  that  x  is  a  function 
of  fand  rand  therefore 

dx  _  dx  dr  _  dx 
dt  dr  dt  4  dr 

This  implies  that  the  slope  or  tangent  of  x(L ,  t )  is 
directly  proportional  to  d.  This  fact  can  then  be  used  to 
adjust  d  to  match  the  mole  fraction  slope. 

This  resulting  identification  procedure  can  be  checked  by 
comparing  oxygen  breakthrough  data  and  additional  nitrogen 
breakthrough  data  to  simulation  predictions. 

In  Table  I,  the  experimental  results  of  nitrogen 
breakthrough  experiments  of  'A,  1,  2,  and  3  atmospheres  are 
presented.  Since  more  than  the  minimum  two  different 
pressures  are  available  (theoretically  all  droop  data  should  be 
identical  at  the  same  pressure),  a  least  squares  fit  of  all  the 
droop  data  was  used  to  find  values  of  *i  and  a.  Note  u*  » 
5/„  where  Afsj  =  28  for  nitrogen  and  M0j  «  32 
for  oxygen.  A  value  of  A,  is  then  found  from  step  2  as 

A ,  =  L(  I  -  <)p0At/  IFofl  -  *i  +  a)tjj 

where  d„  is  the  bed  inlet  gas  density.  From  step  3,  A*  »  (I  - 
n)K ,  Using  slope  adjustment  as  indicated  in  step  4,  the  mass 
transfer  coefficient  d  was  set  at  a  value  of  90  s  " 1  at  1/2  atm 
and  IV',,  -  4.28  gm/s.  This  value  is  adjusted  for  pressure  and 
flow  rate  by  using  a  correlation  due  to  Lighlfoot  [18)  which 
indicates  that  d  is  proportional  to  Wn/P  where  W  is  mass 
flowrate  and  P  is  pressure  In  any  case,  the  value  of  d  did  not 


SB  (NKLE/fraSl 
Fig.  4  Nitrogen  itothonn 


Fig.  5  Oxygon  itothonn 


have  a  large  effect  on  the  characteristic  response.  The  result  of 
this  identification  procedure  gave  the  following  parameter 
identification: 


Ka 

=»  0.153 

*  0.0185 

K , 

-  0.0506 

°k» 

=  0.00611 

b 

-  2.59 

ob 

*  0.31 

are 

the  predicted 

standard 

deviations 

parameters. 

As  a  validation  of  this  procedure  for  identifying  isotherm 
parameters  a  comparison  with  published  isotherms  and 
correlations  is  shown  in  Figs.  4  and  5.  Plotted  in  Fig.  4  are  the 
nitrogen  isotherm  data  of  Peterson  [19),  the  correlation  of 
Ruthven  [10)  based  on  a  statistical  thermodynamics  model, 
and  the  Langmuir  isotherm  identified  in  this  work  in  Fig  5. 
the  oxygen  isotherm  data  of  Wright  [20),  the  correlation  of 
Ruthven  [10),  and  the  linear  isotherm  identified  in  this  work 
are  presented.  The  match  between  the  nitrogen  data  and  the 
model  isotherm  is  excellent.  The  difference  in  the  two  sets  of 
oxygen  data  in  Fig.  5  is  probably  due  to  the  relatively  small 
amount  of  oxygen  absorbed  on  zeolite  5A,  which  requires  a 
very  accurate  gravimetric  measurement.  The  model  identified 
isotherm  bisects  both  sets  of  given  oxygen  data. 

Comparison  of  Simulation  and  Experiment. 

Comparisons  of  the  model,  using  the  identified  parameters, 
with  the  experimental  data  are  shown  in  Fig.  3  and  in  Figs.  6 
through  10.  Shown  in  Fig  3  are  the  data  and  simulation  for 
bed  pressure  at  'A  atm  for  both  oxygen  and  nitrogen  with 
linear  isotherms.  As  can  be  seen  here  for  low  pressure,  a  linear 
isotherm  is  adequate  for  predicting  the  time  responses  The 
same  values  of  A*  and  A*  were  then  used  for  the  3  at 
12 
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Fig  7  '  i  Atmosphere  breakthrough  data  and  simulation  -  nonlinear 

coupled  iaotharma 

mosphere  data,  with  the  results  of  this  simulation  shown  in 
Fig.  6  As  can  he  seen,  the  match  with  experimental  data  is 
\er>  poor  This  indicates  nonlinear  isotherms  are  necessary. 

Comparisons  of  the  model  using  the  nonlinear  nitrogen 
isotherm  with  the  experimental  data  are  presented  in  Figs.  7 
through  10  ,  As  seen  in  these  figures,  the  data  match  has  been 
greatly  improved  over  the  linear  parameter  model. 
Discrepancies  between  model  and  experiment  can  be  at¬ 
tributed  to  experimental  error  and  inaccurately  assumed 
isotherms  The  data  indicate  that  a  nonlinear  oxygen  isotherm 
may  be  appropriate,  but  this  improvement  is  left  for  future 
work 

V  C  (inclusions 

A  J.namic  model  of  the  O.  -  N.  isobanc  breakthrough 
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process  in  a  molecular  sieve  bed  was  developed  and  validated 
experimentally.  This  model  includes  nonlinear  and  coupled 
isotherms.  The  proposed  model  accurately  predicts  the  outlet 
mole  fraction  and  outlet  mass  flow.  In  addition,  isotherms  for 
oxygen  and  nitrogen  were  identified  which  are  in  agreement 
with  previously  obtained  isotherms. 

This  model  is  presently  being  used  in  an  overall  svstem 
model  for  a  pressure  swing  adsorption  process  to  generate 
oxygen  from  aircraft  engine  bleed  air 
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A  Dynamic  Model  of  a  Pressure 
Swing  Oxygen  Generation  System1 

A  full  system  model  for  oxygen  generation  from  a  pressure  swing  adsorption 
process  is  developed.  The  system  under  study  has  two  zeolite  packed  molecular  sieve 
beds  which  are  cycled  with  high  and  low  pressure  in  order  to  produce  oxygen 
enriched  air.  These  On-Board  Oxygen  Generation  Systems  (OBOGS)  are  to  be  used 
for  the  generation  of  aircraft  crew  breathing  requirements.  These  units  are  much 
smaller  in  scale  than  present  industrial  systems  and  are  called  on  to  respond  to 
variable  demand.  These  special  requirements  have  led  to  a  reliable  computer-based 
model  to  simulate  the  relatively  fast  transient  performance  of  OBOGS.  Besides  the 
two  zeolite  beds,  the  system  consists  of  a  rotary  valve  used  to  cycle  the  beds,  purge 
orifice,  and  breathing  plenum.  The  computerized  model  predictions  compare 
favorably  with  experiment.  The  model  is  presently  being  used  in  design  and  control 
synthesis  studies. 


1  Introduction 

As  mentioned  in  (1],  a  promising  system  for  generating 
breathing  oxygen  for  aircraft  crew  uses  a  synthetic  zeolite 
material  to  physically  adsorb  nitrogen  from  engine  bleed  air. 
One  such  system  is  based  on  a  “pressure  swing  adsorption 
cycle"  in  which  nitrogen  is  adsorbed  at  high  pressure  and  then 
desorbed  at  and  exhausted  to  low  pressure.  The  total  process 
produces  oxygen  enrichment  of  air  in  a  cyclical  fashion.  As  a 
first  step  in  the  modeling  of  such  systems,  an  experimentally 
verified  model  of  a  single  bed  system  containing  a  zeolite 
molecular  sieve  bed  was  developed  in  (I).  The  model 
developed  in  |l)  provides  predictions  of  the  bed’s  transient 
outlet  flowrate  and  mole  fraction  for  a  given  inflow  and  mole 
fraction  The  model  was  tested  for  pressures  from  1/2  to  3 
atmospheres  and  inlet  mass  flowrates  of  3.21  gm/s  to  8.56 

gm  s 

Related  work  in  this  area  can  be  found  in  Turnock  and 
Kadlec  |2)  who  study  the  sorption  of  nitrogen  and  methane 
with  pressure  swing  adsorption.  Wright  et  al.  13)  develop  a 
system  model  for  a  thermal  swing  adsorption  cycle  for  CO: 
removal.  Chan  et  al.  (4)  develop  an  equilibrium  theory  for 
pressure  swing  adsorption.  Fernandez  and  Kenney  |5]  study  a 
single  column  pressure  swing  separation  of  oxygen  and 
nitrogen.  Chihara  and  Suzuki  |6]  develop  a  nonisothermal 
pressure  swing  adsorption  model.  Carter  (7)  studies  pressure 
swing  adsorption  for  air  drying. 

Molecular  sieve  zeolite  separation  processes  are  in 
widespread  use  industrially,  but  OBOGS  (acronym  for  On- 
Board  Oxygen  Generation  System)  units  are  much  smaller  in 
scale  and  are  called  on  to  respond  to  a  much  more  variable 
demand  as  determined  by  the  aircraft  crew’s  breathing 
requirements.  These  special  requirements  have  led  to  the  need 
for  a  reliable  computer-based  simulation  tool  that  would 

1  r his  *otk  luppontd  by  Brook!  Air  Foret  B*w.  Aerospace  Medicine.  Crew 

l  echnology 

(  omribuird  hv  ihf  Dvrumic  Sytiemi  and  Control  Division  for  publication  in 
ibe  JoiasAi  o»  Dysahjc  Svstims.  Miajlmmint  and  Control.  Manuscript 
receded  at  ASMfc  Headquarters,  February  19,  1915 


model  the  relatively  rapid  transient  response  in  OBOGS.  This 
model  will  be  useful  for  both  design  modification  and  control 
development. 


2  Description  of  System 

Shown  in  Fig.  1  is  a  schematic  of  a  two  bed  OBOGS 
developed  by  Bcndix  Corporation  for  use  in  military  aircraft. 
The  object  of  the  system  model  is  to  simulate  the  transient 
oxygen  concentration  delivered  to  the  pilot  given  the  pilot’s 
(time  varying)  breathing  mass  flowrate.  This  is  accomplished 
by  mathematically  modeling  the  rotary  valve,  purge  orifice, 
the  two  zeolite  beds,  and  the  breathing  plenum,  and  coupling 
these  separate  models  into  an  overall  system  model. 

The  pressure  swing  adsorption  cyclic  operation  of  the  beds 
can  be  described  as: 

1.  The  inlet  of  bed  I  is  supplied  air  at  supply  pressure 
through  the  rotary  valve,  while  bed  2  is  being  purged  to 
exhaust  pressure  through  the  valve.  During  this  stage 
breathing  flow  is  obtained  from  bed  I . 

2.  As  the  rotary  valve  turns,  the  inlet  of  bed  1  is  switched  to 
exhaust  pressure  and  the  inlet  of  bed  2  is  switched  to  supply 
pressure.  During  an  initial  transition  stage,  the  pressure  in  bed 
1,/V  is  still  higher  than  the  pressure  in  bed  2 ,P:.  In  this  case, 
there  is  flow  out  of  both  sides  of  bed  I  and  into  both  sides  of 
bed  2.  Breathing  flow  will  come  from  the  higher  pressure  bed 
I. 

3.  As  the  pressure  drops  in  bed  I  due  to  outflow  and  builds 
in  bed  2  due  to  inflow,  P ,  will  become  greater  than  P , .  This  is 
the  identical  case  of  (I)  above,  but  with  bed  2  the  higher 
pressure  side.  The  stages  I,  2,  and  3  above  describe  1/2  of  a 
production  cycle.  The  other  half  cycle  is  symmemc  to  the  first 
half  cycle. 

The  components  of  the  system  are  the  zeolite  molecular 
sieve  beds,  rotary  valve,  bypass  or  purge  valve  between  the 
<wo  beds,  and  the  breathing  plenum.  These  components  are 
modeled  as  discussed  below. 
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changes  in  the  free  stream  concentration,  and  isobartc  beds 
(the  pressure  drop  across  the  beds  is  small),  the 
oxygen/nitrogen  adsorption  in  the  beds  can  be  modeled  as 


3(1*4) 

3? 


=  A0(y/C-  I) 


(1) 


3  (yA) 
3z 


=AB(yA  -y/u) 


(2) 


Fig.  2  Concentric  Molite  bed 


2. 1  Zeolite  Bed  Model.  Concentric  beds  were  used  in  the 
experimental  study  and  one  is  depicted  in  Fig.  2.  As  discussed 
in  (I |,  if  we  assume  one  dimensional  flow,  negligible  axial 
diffusion,  diffusion  into  the  zeolite  described  with  a  mass 
transfer  coefficient,  convective  terms  dominate  temporal 


^--dCr-O  0) 


=  -d(yA  -CA)  (4) 

where 

A  *  bed  cross-sectional  area 

overall  mass  transfer  coefficient 
£  =  </(  1  (-void  f  ran  ion  of  gas  volume  to  total 

volume 

C  =  CA  +  CB  total  gas  phase  concentration 
CA  =  oxygen  gas  phase  concentration 
y=yA+yt  total  interfacial  concentration 
yA  =*  oxygen  interfacial  concentration 
u  *  convective  velocity 

rr*nA  +  nB  total  adsorbed  phase  concentration 
nA  *  oxygen  adsorbed  phase  concentration 
y  =  uCA  convective  oxygen  flux 

Due  to  the  relatively  fast  heat  transfer,  it  was  observed  that 


Nomenclature 


A  =  bed  cross-sectional  area,  mJ 
Ad  =  orifice  area,  mJ 
A,  =  inner  annulus  area,  m 2 
A„  =  outer  annulus  area,  m] 

C  =  gas  concentration,  kgmole/m 5 
C„  =  valve  discharge  coefficient 
C e  =  breathing  plenum  gas  concentration, 
kgmole/m1 

d  =  total  mass  transfer  coefficient.  I /s 
k  =  ratio  of  specific  heats 
L  =  bed  length,  m 

M  =  average  molecular  weight  of  gas  in  plenum, 
kg/kgmole 

n  =  adsorbed  phase  concentration,  kgmole/m} 
PU.P,,P:.P,  =  supply,  Bird  I,  Bed  2,  exhaust  pressures, 
(N/nv) 

Q  -  volume  flowrate,  mJ/s 
i  -  time,  s 

u  =  convective  velocity,  m/s 
H  =  mass  flowrate,  kg/s 
y  =  molecular  flux  of  oxygen,  kgmole/m;-s 


z  =  bed  distance,  m 
0  =  d(\  -*)/«  (1/s) 

y  =  interfacial  concentration,  kgmole/m* 
A/  =  time  step,  s 
Az  =  space  lump,  m 
«  *  void  fraction 
A  =  denotes  property  of  oxygen 
B  *  denotes  property  of  nitrogen 
BR  *  denotes  breathing  property 
i  =  space  index 
j  *  time  index 
L  *  denotes  property  at  z  =  1 
P  =  denotes  purge  orifice  property 
0  -  denotes  property  at  z  =  0 

1  »  denotes  bed  I  property 

2  -  denotes  bed  2  property 

Superscripts 

-  *  denotes  property  just  before  z  =  L/2 
*  -  denotes  property  just  after  z  =  L/2 
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experimental  beds  remained  at  essentially  constant 
ernperature  during  a  cycle  even  though  the  adsorption 
process  is  not  intrinsically  isothermal.  A  more  detailed  study 
ot  the  thermal  properties  of  the  system  is  the  subject  of 
ongoing  research.  In  the  isothermal  case,  >4  and  yt,  the 
imcrfacial  concentrations  of  oxygen  and  nitrogen,  respec¬ 
tively.  are  related  to  nA  and  nB,  the  adsorbed  oxygen  and 
nitrogen  concentrations,  by  the  equilibrium  adsorption 
isotherms  |8)  As  discussed  in  {1].  these  isotherms  can  be 
expressed  as 

>4=>4('ix.'t»)  (5  a) 

•>a  =  ya<««."»)  (5b) 

tor  nonlinear  coupled  isotherms.  The  pure  component 

isotherms  were  identified  as 

>'<=K4n*4  (6fl) 

-n\/b)  (6b) 


»uh  A,  =0.153  kg  mole  gas  phase/kgmole  adsorbed  phase, 
A  „  =0  0506  kgmole  gas  phase/kgmole  adsorbed  phase,  and 
b  -  2.59  kgmole  adsorbed  phase.  The  coupling  is  predicted  by 
using  an  ideal  adsorption  method  (I). 

The  diffusion  coefficient  d  was  empirically  chosen  as  200 
I  s  in  order  to  best  match  the  data.  This  value  is  higher  than 
tne  values  identified  in  (1).  This  is  probably  due  to  the  high 
velocities  in  the  beds  during  exhaust  and  initial  intake  of  air. 

2.2  Valve  Model.  The  orifices  of  the  rotary  valve  and 
the  purge  orifice  are  modeled  as  standard  compressible  flow 
pneumatic  orifices  as  described  in  (9).  The  basic  equation  for 
the  flow  of  a  gas  through  an  orifice  is  of  the  following  form 

(7) 

where 

*•  =  mass  flowrate 
.4,,  is  orifice  area 

C.  is  a  dimensionless  discharge  coefficient 
f.  is  upstream  stagnation  temperature 
is  upstream  stagnation  pressure 
P ,  is  downstream  pressure 
/ is  a  function  of  P„.  T„,  and  P^/P*. 

for  pressure  ratios  Pd/Pu  greater  than  the  choked  flow  ratio 
and  less  than  1 .  the  funciion/ can  be  expressed  as 

J=C.  ~~  (PJPJ'  ‘(I -(PjP^t  "‘)' J;  (8) 

for  pressure  ratios  less  than  the  choked  flow  ratio. 


where  C,  and  C:  are  functions  of  k,  the  ratio  of  specific 
heats,  and  the  gas  constant.  These  parameters  are  taken  as 
constants  for  this  study.  For  the  rotary  valve  the  orifice  area* 
are  cyclic  functions  of  the  rotation  angle  of  the  valve.  The 
discharge  coefficients  were  chosen  between  values  of  0.6  and 
1  in  order  to  best  match  the  measured  experimental  flow. 

2.3  The  Breathing  Plenum.  The  breathing  plenum  is 
modeled  as  an  ideal  isobaric  two  stage  mixing  chamber.  The 
transfer  function  between  inlet  mole  fraction  and  outlet  mole 
traction  can  be  expressed  as: 


=  1 

x,  («+  t)J 

where 


(9) 
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Fig.  1  Conceptual  bond  graph  ot  ayatent 


is  breathing  mass  flowrate.  Cp  is  the  plenum  molecular 
concentration,  VM  is  approximately  half  of  the  plenum 
volume.  A/  is  an  average  molecular  weight.  In  dertvmg  these 
equations  the  following  assumptions  have  been  made: 

1 .  The  plenum  remains  at  approximately  constant  pressure 
and  therefore  constant  Cr.  This  is  not  strictly  true  as 
breathing  oxygen  is  being  switched  between  the  bedc.  but  for 
our  prototype  system  this  assumption  is  reasonable. 

2.  The  molecular  flux  into  and  out  of  the  plenum  can  be 
found  approximately  from  the  mass  flowrate  and  an  average 
molecular  weight.  This  assumption  is  quite  reasonable  for 
oxygen  and  nitrogen  since  their  individual  molecular  weights 
do  not  differ  substantially. 

3.  The  mixing  volumes  for  the  two  idealized  chambers  are 
equal  and  can  be  found  from  half  the  plenum  volume  plus  a 
correction  for  mixing  in  the  piping. 

3  Solution  Technique 

3.1  Causal  Considerations.  The  overall  system  model 
can  be  depicted  by  the  word  bond  graph  (I0|  shown  in  Fig.  3. 
Besides  the  physical  interconnections,  the  graph  also  indicates 
the  important  causal  relationships  in  the  system  model.  The 
causality  shown  in  the  graph  can  be  read  as 

1.  Pressures  P0  and  P,  are  supplied  as  inputs  to  the  rotary 
valve  while  the  valve  determines  the  flowrate  Q0I  from  these 
two  pressures. 

2.  P !  is  determined  from  bed  I  which  has  inputs  the  flow 
Go.  *nd  the  flow  QLI . 

3.  Pressures  P ,  and  Pi  are  the  inputs  to  the  purge  valve  and 
the  valve  supplies  the  flowrates  Gn  and  QK  to  the  system 
Qtll,  and  G*«  are  determined  by  the  breathing  flow  and  bed 
pressures  in  the  check  valve/  plenum  system  as  described  in  the 
next  section. 

4.  P:  is  determined  by  bed  2  which  has  inputs  the  flows  Go: 
and  Gci- 

5  Giz  **  determined  by  the  rotary  valve  with  inputs  to  the 
valve  of  P:  and  P, . 

In  addition  to  these  causal  statements,  the  spatial  molecular 
concentration  of  the  two  gases  in  the  beds  is  determined  in  the 
direction  of  the  flowrate  i.c.,  the  concentration  must  be 
determined  from  upstream  conditions.  This  will  be  discussed 
more  thoroughly  in  the  next  section. 

The  causality  described  above  and  shown  on  the  graph  is 
integral  causality  for  the  beds  if  they  are  treated  as  capacitive 
fields  |ll]  or  mass  storage  devices.  Using  other  causal 
arrangements  is  possible  but  will  result  in  complicated  and 
computationally  time  consuming  iteration. 

3.2  Solution  Algorithm.  Given  input  supply  pressures 
P0  and  Pi,  breathing  flow  and  initial  adsorbed  phase 
concentrations  rt U.0)  and  n,  (:,0)  in  the  two  beds,  the  system 
equations  can  be  solved  as  a  boundary  value  problem  in  space 
and  a  marching  problem  in  time  The  solution  technique 
consists  of: 

(0  Given  ir(z.O)  and  »4U.O),  the  tnterfactal  gas  phase 
concentrations  >U,0)  and  >  4  (z,0)  can  be  calculated  Irom  the 
equilibrium  isotherms  (5).  Equation  ( I )  can  then  be  integrated 
along  the  enure  length  of  each  bed  as 

<Ci  -  Go  ,3Ls2  +  A„JL  2)C 

r  i  :  r 1 

yd:  +  A3\  ydz  <101 

Jo  J  i  i 
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Fig.  4  0»yg*fi  mol*  If  action  *-w»u»tim*  (10  (liumtfi)bf»4thingtk>w)  Flg.S  Oxygen  mole  traction  versus  IMw  (SO  (lil/min)  breathing  flow) 


Fig  s  OinnuiwMInettMiMuiUiNlN (KIMfi) bnMMif  On)  Fig.  7  Oxygen  mole  traction  tmn tune (70 (IttfiMn)  breelhing  Now) 


lor  bed  I.  Qn  =  Qo,{P„.P,)  where  P,  is  proportional  to  C, , 
the  gas  phase  concentration  in  bed  1 ,  for  isothermal  beds  and 
assuming  a  perfect  gas.  and  Qt  =Qt  (/*, .  P2)  -  QM,  where 
P:  is  proportional  to  C<,  the  gas  concentration  in  bed  2  and 
Qt, i  accounts  for  the  volume  flow  out  of  bed  1  due  to  crew 
breathing  For  product  flow  out  of  bed  I.  Qlfl  is  propor¬ 
tional  to  for  an  isobaric  and  isothermal  phenum.  For 
bed  2.  Q{  -Qli:[P,.  P.)  -  Qim:  and  QL  =  Q^iP^P,) 
where  QMlt:  accounts  for  breathing  flow  out  of  bed  2.  Note 
that  either  or  Qltl  will  be  zero  depending  on  which  bed 
had  the  lower  pressure.  Equation  (10)  written  out  for  both 
beds  represents  two  coupled  nonlinear  algebraic  equations  for 
the  i*o  unknowns  C,  and  C2.  These  equations  can  be  solved 
by  an  iteration  procedure  consisting  of  (I)  assuming  con¬ 
centration  C:  and  solving  the  resultant  single  nonlinear 
algebraic  equation  for  C,  using  a  modified  false  position 
algorithm.  <2)  using  this  value  of  C,  to  calculate  a  new  value 
of  C  ,  (i)  iterate  until  the  solution  procedure  converges. 
Consergence  of  this  procedure  requires  on  the  order  of  three 
to  four  iterations  The  integrations  involved  in  (10)  are 
performed  numerically  by  a  trapezoidal  rule,  and  are  only 
done  once  per  time  step.  This  iteration  procedure  amounts  to 
a  solution  of  the  boundary  value  problem  in  the  system. 
in)  Spatially  integrate  equation  <1 )  numerically  as 

u  ~u  ,  +d<>, +y  _  ,)/2C-dAz,  (II) 

where  i  is  the  index  (or  the  space  dimension,  c  «=  rAz.  and  Az  is 
the  spjse  lump  sire  This  integration  is  done  in  two  steps  from 
0  to  l  2  and  then  from  1/2  to  l  with  condition  uU'  )  » 
It;  uric  )/-!(;*)  at  ;  =  £'  2  to  account  for  the  area 


change  in  the  outer  annulus  and  inner  cylinder  of  the  bed 
Note  that  the  starting  velocity  u0*Q0/4( 0)  is  known  from 
step  (f). 

(iii)  Using  the  velocity  distribution  calculated  in  (ii). 
equation  (2)  is  a  first  order  linear  ordinary  differential 
equation  in  space  for  fixed  time.  Approximate  the  velocity  u 
and  the  interfacial  concentration  >  ,  over  a  space  lump  as 

u*r0  +  r,z 

( r  -  I  lAz  <  z  <  /Az 
y, 
where 

r0  »  u,  r,  =  <ir  -u,  ,)/Az 
m  Ta, - 1 .  r,  «(*4,  -y«,  ,)/Az 

Using  these  approximations,  equation  (2)  can  be  iniegr ated  m 
closed  form  as  shown  in  |l|  The  direction  of  the  integration 
will  always  be  in  the  direction  of  the  velocity  flow  Which,  as 
described  in  Section  2.  might  be  (I)  flow  in  one  side  of  the  bed 
and  out  the  other.  (2)  flow  into  both  sides  of  the  bed.  or  (3) 
flow  out  of  both  sides  of  fhe  bed  In  cases  (2)  and  (3).  there 
will  be  at  least  one  point  in  the  bed  that  has  zero  flow  As  the 
velocity  in  the  bed  is  found  as  described  in  (ii)  above,  the  point 
of  zero  velocity  is  stored  For  case  (I),  the  integration  will 
proceed  in  the  direction  of  flow,  lor  case  (2)  the  integration 
will  proceed  from  the  sides  of  the  bed  to  the  point  of  zero 
flow;  and  for  case  (3)  the  integration  will  start  at  the  point  of 
zero  velocity  and  proceed  out  the  two  sides  Although 
physically  possible,  multiple  points  of  zero  flow  will  nol  lKkUr 
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Fig  (  Oxygen  mol*  traction  versus  cycle  tlmo  (30  (llt/mtn)  breathing 
•low) 

in  normal  operation  of  OBOGS  and  this  case  is  not  included 
ai  present  in  ihe  computerized  solution. 

The  ability  to  solve  equation  (2)  in  closed  form  is  an  im¬ 
portant  numerical  result.  The  oxygen  profiles  in  the  beds  can 
be  very  steep,  and  if  a  completely  numerical  procedure  is 
applied  to  (2),  it  will  require  small  increments  for  Az  in  order 
to  be  stable  and  accurate.  This  will  result  in  very  long  com¬ 
putation  times  even  on  a  high  speed  computer . 

(it  i  Equations  (3)  and  (4)  are  then  integrated  in  time  via  a 
simple  Euler  method  as 

n,  ,  ;.| -rf(r -OAZ  (12) 

, -d<7a, -Ct,)Af  (U) 

where  j  is  the  time  index.  and  A/  is  the  time  step. 

(t  I  The  input  oxygen  mole  fraction  x ,  to  the  breathing 
plenum  is  found  from  bed  I  at  z  =  L  or  bed  2  at  z»0 
depending  on  which  bed  has  higher  pressure.  The  mole 
fraction  from  the  higher  pressure  bed  is  computed  as  CA/C. 
This  input  mole  fraction  is  used  with  equation  (9)  and  the 
breathing  (low  W  to  solve  for  the  output  mole  fraction,  x0. 
Equation  (9)  is  solved  as  two  coupled  ordinary  differential 
equations  via  a  simple  Euler  method. 

(til  March  in  time  by  proceeding  to  step  </>  for  the  new 
time  level 

The  time  step  and  space  lumps  were  chosen  by  performing 
numerical  experiments  to  insure  accuracy  and  stability  of  the 
solution  procedure.  The  model  simulation  time  is  ap¬ 
proximately  twice  that  of  real  time  when  computed  on  a  CDC 
t  sber  '50  1 10 

4  C  omparison  of  Model  With  Experiment 

The  prototype  experimental  system  is  depicted  in  Figs.  1 
and  2  The  concentric  moleculat  sieve  beds  have  an  inner  area 
of  3ft  5  cm  and  an  effective  outer  area  of  136.  cm2.  The 
effeeme  flow  length  of  a  bed  is  39.4  cm.  The  beds  are  loaded 
with  l  mde  *  A  sorbent  with  zeolnte  panicles  of  approximately 
140  io  «33  microns  in  diameter  The  void  fraction  of  the  beds 
"  i  -  0  1  All  experiments  were  run  at  room  temperature  of 
’’’C  No  noticeable  variation  in  temperature  was  observed 
during  ihe  system  operation  which  implies  the  isothermal 
assumption  is  valid  for  the  system  under  study  The  fully  open 
supply  pressure  valve  has  an  effective  diameter  of  .78  cm.  and 
a  discharge  coefficient  of  0  8  The  fully  open  exhaust  pressure 
salve  has  an  effective  diameter  of  I  12  cm  and  a  discharge 
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coefficient  of  0.8.  The  purge  orifice  has  a  diameter  of  0  191 
cm  and  a  discharge  coefficient  of  0.75.  The  breathing  plenum 
volume  it  .0018  m’  an  additional  volume  of  .0006  m1  is  added 
to  this  quantity  to  account  for  mixing  in  the  connections. 

In  this  paper,  comparison  with  of  the  model  is  made  with 
three  basic  types  of  experiments:  (I)  transient  startup.  (2) 
cycle  time  variation,  and  (3)  change  of  active  volume  of 
zeolite  in  the  beds. 

4.1  TraasSeat  Startup  ol  OBOGS.  Shown  in  Figs  4 
through  7  are  the  results  from  the  experimental  data  and  the 
simulation  model  for  the  OBOGS  at  initially  atmospheric 
pressure,  the  rotary  valves  turned  on  at  time  *  0.  a  supply 
pressure  of  0.31  Mpa.  an  exhaust  pressure  of  .10  Mpa.  and 
cycle  time  of  10.7  seconds  These  data  are  shown  for  four 
different  constant  breathing  flowrates  of  0.214,  0  642,  t  07, 
and  1 .50  gm/rain  (10,  30,  50,  and  70  standard  liters,  mm)  As 
can  be  seen  in  these  figures,  the  model  does  an  excellent  job  of 
matching  the  transient  experimental  data  for  the  breathing 
oxygen  mole  fraction  over  the  range  of  breathing  flowrates 
The  increased  ripple  in  the  data  and  simulation  for  increased 
flow  is  due  to  less  time  for  complete  mixing  in  the  plenum  as 
the  breathing  flow  goes  up.  The  model  describes  this  effect 
very  well. 

These  plots  indicate  a  method  for  controlling  the  outlet 
mole  fraction.  This  is  desirable  since  both  high  and  low 
oxygen  mole  fraction  can  cause  physiological  problems  (121  to 
the  aircraft  crew.  By  dumping  product  gas  to  artificially 
increase  the  breathing  flow,  the  outlet  mole  fraction  can  be 
lowered  if  too  high.  If  the  mole  fraction  is  too  low.  product 
gas  dumping  can  be  decreased. 

4.2  Effect  of  Cyrie  Time  Variation.  In  Fig.  8,  the  ex¬ 
perimental  and  simulation  results  for  a  variation  in  cycle  time 
are  shown.  As  can  be  seen,  the  model  does  a  good  )ob  of 
describing  changes  in  oxygen  output  wuh  cycle  lime  variation 
The  results  in  this  figure  also  demonstrate  a  method  lor 
controlling  the  outlet  oxygen.  By  simply  varying  the  cvcle  time 
of  the  process,  control  of  oxygen  output  is  possible  T  he  c  v  Je 
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rime  variation  can  be  easily  implemented  by  using  a  variable 
vpejd  drive  for  the  rotary  valve.  This  control  scheme  has  the 
'  advantage  of  not  using  any  extra  product  gas  which  might  be 
■  desirable  for  the  contaminant  filtering  properties  of  the  beds. 
‘Tha  study  of  this  control  scheme  is  the  subject  of  ongoing 
research  and  these  results  will  be  presented  at  a  later  date. 

4.3  Bed  Volume  Variation.  Figure  9  shows  the  effect  of 
changing  the  volume  of  \he  beds,  this  was  done*  ex¬ 
perimentally  (13)  by  filling  the  Beds  with  fractional  amounts 
of  molecular  sieve  which  had  been  deactivated  with  water. 
From  simulation  results,  itVas  found  that  only  bed  volume 
and  not  shape  is  of  importance  for  oxygen  output.  This  might 
be  an  important  design  result,  fo(  it  implies  that  the  shapes  of 
the  beds  are  not  critical.  This  means  that  the  beds  could  be 
shaped  for  space  considerations  in  the  aircraft. 

Another  interesting  result,  is  the  relative  insensitivity  of  the 
beds  to  volume  changes  until  a  critical  volume  is  reached 
(around  75  percent  for  these  beds).  This  also  has  important 
design  implications.  It  means  that  siring  the  system  for  in¬ 
creased  flow  is  not  just  a  simple  increase  in  bed  size,  but  a 
combination  of  volume  changes  and  other  design  parameters. 
In  particular,  the  valve  areas  and  discharge  coefficients  are 
important  design  values. 

In  addition,  these  results  indicate  a  potential  problem  in 
operation  If  the  beds  were  to  become  contaminated  with 
water,  which  has  the  effect  of  decreasing  the  effective  volume 
ol  the  beds,  there  would  be  no  outward  indication  of  per¬ 
formance  degradation  until  a  critical  volume  were  reached, 
and  at  that  point  the  beds  would  rapidly  lose  performance 
wuh  any  further  contaminarion. 

5  Conclusions 

A  relatively  simple  model  of  a  pressure  swing  oxygen 
generation  system  was  developed  in  this  study.  The  com- 
phtenred  model  includes  the  zeolite  beds,  valves,  and  plenum 
breathing  chamber  in  the  system.  The  model  predictions 
compare  favorably  with  detailed  experimental  data.  Ongoing 


work  in  this  area  include  the  prediction  of  temperature  and 
contaminant  effects  and  control  of  output  breathing  oxygen. 
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COMPUTER  MODEL  FLOWCHART 

I  (inventing  Equations 


(I) 


fXUC) 

()/ 


|)<Y-C> 


p-a-('-a 

£ 


(2) 


U(UCA) 

()/. 


P<?A  -  cA) 


(3) 


''nA 

—  -d(yA-cA) 

<n  A  A 


(4)  ^  -d(y-C) 

dt 


II  Procedures 

i  1 1  initial  guess  P(  I )  and  P( 2 ) 
i  2)  l  sc  adsorption  isotherms  to  calculate  yA  and  yB 


(  '  i  Integrate  equation  ( 1 )  to  get  C: 


(I  A  r  A  )C 

in  t)  I 


(y-C)dz 


( 


(l  ,  A 


—  Cldz 

«>  I 


J.\ 
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values  will  be  consistent  with  the  supply  pressure  and  exhaust  pressure  at  each 
time  step. 

(5)  Substitute  P(l)  and  P(2)  (calculated  in  step  4)  into  the  valve  eq.  can 
get  w(l),  w(2) ,  and  w(3) 

(6)  Since  Ui  «  W  /  (Ci*Ai*Mi) ,  where  Mi  is  the  molecular  weight,  U0,  ,U||  ,  and 
U02 and  can  be  calculated. 

(7)  Assume  velocity  is  linear  along  the  bed  length.  By  integrating  eq.  (2), 
CA(Z)  can  be  obtained. 

(8)  Use  eq.  (3)  and  eq.  (4)  to  get  new  values  for  ^ , %  and  7  for  each  time 
step. 

(9)  T  =  T  +  dT. 

Get  a  new  valve  area  and  repeat  the  procedure  again. 


III.  Flowchart  of  the  computer  program 


(1)  Sign  Convention 


(2)  Flowchart 


Please  see  attached. 
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update  time  to  calculate  new  and  n 


sets  initial  values  for  P(I),P(2),  w(l),w(2) 
w(3) ,x(l) ,x(2) ,x(3) ,  nl  and  n2 


solves  for  w(l) ,  w(2),w(3),  P(1),P(2), 

U6i ,  Uoi,  U02  ,  and  U»2 

where  U<a  and  Un  are  the  velocities  at  the 
end  of  bed  #1.  U<o  and  Uf2  are  referred  to 
bed  n. 


''T* - —1  CALL  CONCA 

i  [ - 


solves  for  velocity  distribution,  l'(Z,t), 
and  mole  fracrion  CA(Z,t)  of  each  bed 


on-line  ReGIS 
plotting  for 
XBR,  and  w's 


APPENDIX  C 


25 


APPENDIX  C 


The  inDuts  to  the  orogram  are,  as  in  the  menu: 


SYSTEM  PARAMETERS 

1.  supply  pressure 

at 

40.00 

PSIA 

2.  outlet  pressure 

= 

15.00 

PSIA 

3.  final  observation  time 

= 

30.00 

SEC. 

4.  cycle  time 

as 

10.70 

SEC. 

5.  breathing  flowrate 

s 

10.00 

STD  LIT/MIN 

GEOMETRIC  PARAMETERS 

6.  by-oass  valve  diameter 

= 

0.0750 

IN. 

7.  supply  valve  diameter 

m 

0.3062 

IN. 

8.  outlet  valve  diameter 

m 

0.4387 

IN. 

9.  bed  length 

K 

15.50 

IN. 

10.  outer  bed  diameter 

m 

5.7300 

IN. 

11.  inner  bed  diameter 

a 

2.1800 

IN. 

BED  PARAMETERS 

12.  KA 

= 

0.1530 

KGMOLE  02  ABSORBED/KGMOLE  GAS 

13.  KB 

« 

0.0490 

KGMOLE  N2  ABSORBED/KGMOLE  GAS 

14.  B 

m 

2.6500 

15.  diffusion  coefficient 

s 

200.00 

1/SEC 

16.  void  fraction 

- 

0.37 

The  user  has  the  following  oDtions  after  he  has  reviewed  the  menu: 

1.  change  parameters 

2.  insert  a  step  change  in  the  breathing  flowrate,  or 

3.  exit  the  program 


First,  changing  the  parameters  is  done  by  entering  the  corresponding  parameter 
number  appeared  in  the  menu.  The  program  then  prompts,  for  examnle: 

ENTER  THE  PARAMETER  AND  PRESS  "RETURN" 


SUPPLY  PRESSURE  (PSIA) 


CURRENT  VALUE 
40.00 


Once  the  new  value  has  been  assigned,  the  program  goes  back  to  the  menu  with 
modified  parameter(s)  and  waits  for  the  next  input. 


If,  on  the  other  hand,  the  user  wishes  to  insert  a  step  change  in  breathing 
flowrate,  he  enters  the  number  for  that  option  (which  is  "59"  in  this  program). 
The  program  will  ask  for  the  starting  and  ending  time,  as  well  as  the  value 
for  the  step  change  and  again  goes  back  to  the  menu. 


After  the  completion  of  modifications  on  the  parameters,  the  user  enters  "0” 
to  run  simulation.  The  program  responds  with: 
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THE  FOLLOWING  PLOTS  ARE  AVAILABLE  TO  VOU  ON  THE  TERMINAL 


1  MOLE  FRACTION  OF  OXVCEN  VS. TIME  (DATA  FILE=MOLEFRAC  DAT) 

2  INLET  MASS  FLOWRATE  VS  TIME  (DATA  F I LE=FLOWR ATE  DAT > 

3  OUTLET  MASS  FLOWRATE  VS  TIME  (DATA  F ILE=FLOWRATE  DAT) 

4  A  DYNAMIC  SIMULATION  FOR  ONE  CYCLE  (OXVCEN  VS.  TIME, 

DATA  FILE=PROFILE  DAT) 


PRESS  THE  CORRESPONDING  «  TO  HAVE  THE  OUTPUT  SHOWN  AS  THE 
SIMULATION  TRANSPIRES 


I  Ik-  ’.  >ei  chooses  the  desired  output,  and  the  program  gives  the  following: 


1  ENTER  "0"  TO  RUN  SIMULATION 

2  ENTER  "IS"  TO  CHANGE  THE  PARAMETERS 

3  ENTER  "99"  TO  EXIT  PROGRAM 


LMER  THE  CORRESPONDING  #  AND  PRESS  "RETURN" 


•  1 1 •* I  veils  lor  i he  decision  from  the  user.  That  is.  it  is  still  possible  for 
i in  usei  to  go  back  to  the  menu  before  he  nins  the  simulation. 

hr.  ills .  the  program  asks  for  a  last  item: 


1C'  'OU  WISH  10  DOUBLE  THE  NUMBER  OF  SPACE  LUMPS  IN 
H  -EMULATION  FOR  BETTER  ACCURACY  ■>  (  CURRENT  NUMBER  USED 

,c.  SO  )  (V/N) 


I  .in  Hue  to  the  t.ni  ili.it  when  using  one  hundred  lumps  the 
■  i  ii  .it i. mi  is  time  consuming  for  the  low  breathing  flowrate  cases 
I  •  Uetauli  is  thus  set  to  be  M)  tor  faster  simulation 


MS 


CMS 
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:  URFEMT  SVSTEM  PARAMETERS 

ARE: 

1 

SUPPLY  PRESSURE= 

40.00 

PSIA 

OUTLET  PRESSURE= 

15.00 

PSIA 

FINAL  OBSERVATION  TIME= 

30.00 

SEC 

, 

CVCI.E  TIME= 

10.70 

SEC 

3PEATHING  FL0WRATE= 

10.00 

STD  LIT/MIN 

TORRENT  GEOMETRIC  PARAMETERS  ARE: 

3V-PASS  VALVE  DIAMETER  = 

0.0750 

IN 

SUPPLY  VALVE  DIAMETER  = 

0.3062 

IN 

OUTLET  VALVE  DIAMETER  = 

0.4387 

IN 

■j 

BED  LENGTH  = 

15.5000 

IN 

1-:! 

OUTER  BED  DIAMETER  = 

5 . 7300 

IN 

1  1 

INNER  BED  DIAMETER  = 

2.1800 

IN 

THE  CURRENT  BED  PARAMETERS 

ARE: 

1 

KA  •: 

0.1530 

KGMOLES 

02  ABSORBED /KGMOLE 

CAS 

1  i 

KB  - 

0 . 0490 

KGMOLES 

N2  ABSORBED /KGMOLE 

GAS 

14 

B  = 

2 . 6500 

15 

DIFFUSION  COEFFICIENT  = 

200 . 00 

1/SEC 

16 

VOID  FRACTION  = 

0.370 

ENTER  CORRESPONDING  PARAMETER  #  TO  CHANGE  PARAMETER 
ENTER  "O'1  TO  SIGNIFV  PARAMETERS  ARE  CORRECT — RUN  SIMULATION 
ENTER  ""59"  TO  INSERT  A  STEP  CHANGE  IN  BREATHING  FLOWRATE 
ENTER  "09"  TO  EXIT  PROGRAM 
1 


MM  MM  I  MMMMMMMMM************************* 

MEN  THE  PARAMETER  AND  PRESS  "RETURN" 


CURRENT  VALUE 


PRESSURE  (PSIA) 


40  0000 


28 


ft 


. .  , 
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CURRENT  SYSTEM  PARAMETERS  ARE 


1 

SUPPLY  PRESSURE3 

30  00 

PSIA 

OUTLET  PRESSURE3 

15  00 

PSIA 

-INAL  OBSERVATION  TIME3 

30  00 

SEC 

", 

CYCLE  TIME3 

10  70 

SEC 

BPEATHINC  FLOWRATE3 

10  00 

STD  LIT/MIN 

CURRENT  CEOMETRIC  PARAMETERS  ARE 

BY-PASS  VALVE  DIAMETER  3 

0  0750 

IN 

SUPPLY  VALVE  DIAMETER  3 

0  3062 

IN 

OUTLET  VALVE  DIAMETER  3 

0  4387 

IN 

3ED  LENCTH  3 

15  5000 

IN 

]-■ 

OUTER  BED  DIAMETER  3 

5  7300 

IN 

1  ; 

INNER  BED  DIAMETER  3 

2  1800 

IN 

) HE  CURRENT  BED  PARAMETERS 

ARE: 

t  _ 

KA  ■= 

0  1530 

KGMOLES 

02  ABSORBED/KGMOLE 

GAS 

1  • 

KB  == 

0  0490 

KGMOLES 

N2  ABSORBED/KGMOLE 

GAS 

1". 

B  = 

2  6500 

1 

DIFFUSION  COEFFICIENT  3 

200  00 

1/SEC 

1‘. 

VOID  FRACTION  3 

0  370 

ENTER  CORRESPONDING  PARAMETER  *  TO  CHANCE  PARAMETER 
ENTER  "O'  TO  SIGNIFY  PARAMETERS  ARE  CORRECT — RUN  SIMULATION 
ENTER  "59"  TO  INSERT  A  STEP  CHANGE  IN  BREATHING  FLOWRATE 
ENTER  "99"  TO  EXIT  PROCRAM 

59 


T  t  ltt.tttttttt******t*t***t**t******ttt***t**t**** 

ENTER  THE  TIME  AT  WHICH  THE  STEP  CHANCE 
IN  ERE A  THING  FLOWRATE  IS  TO  OCCUR 

THE  FINAL  OBSERVATION  TIME  IS  CURRENTLY  30.00  SECONDS 


1 


! 


.  'tut'  U*itt*ii:****tt*ttt****t*t**t****tt****t** 

",  ■>  THE  TIME  AT  WHICH  THE  STEP  CHANCE 
)  CLAIMING  FLOWRATE  IS  TO  END 

•  ! N'ti.  OBSERVATION  TIME  IS  CURRENTLY  30  00  SECONDS 
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t ,*tit**t**t*tt*******tt**t*******tt*****t*ttt*tt 

ENTER  THE  NEW  BREATHING  FLOWRATE  OCCURING  AS 
A  STEP  CHANGE  FROM  T=  10.00  TO  T=  20.00  SECONDS 

THE  CURRENT  BREATHING  FLOWRATE  IS  10.00  STD  LIT/MIN 


:?(■> 


!  i n*tit*****************t*tt**t*****t*****t*t***t 


I  HE  BREATHING  FLOWRATE  OF  30.00  STD  LIT/MIN 
WILL  BE  INPUT  AS  A  STEP  AT  T=  10.00  SECONDS 
'  NI)  WILL.  END  AT  T=  20.00  SECONDS 


IJRPENT  SVSTEM  PARAMETERS 

ARE: 

1 

SUPPLV  PRESSURE3 

30  00 

PSIA 

> 

OUTLET  PRESSURE3 

15.00 

PSIA 

• 

FINAL  OBSERVATION  TIME3 

30.00 

SEC 

t 

CVCLE  TIME3 

10  70 

SEC 

BREATHING  FLOWRATE3 

10  00 

STD  LIT/MIN 

I  URREN1  GEOMETRIC  PARAMETERS  ARE: 

6 

BV-PASS  VALVE  DIAMETER  3 

0  0750 

IN 

7 

SUPPLV  VALVE  DIAMETER  3 

0 . 3062 

IN 

8 

OUTLET  VALVE  DIAMETER  3 

0  4387 

IN 

9 

BED  LENGTH  3 

15.5000 

IN 

10 

OUTER  BED  DIAMETER  3 

5 . 7300 

IN 

1  1 

INNER  BED  DIAMETER  3 

2  1800 

IN 

'HE  CURRENT  BED  PARAMETERS 

ARE: 

12 

KA  = 

0  1530 

KGMOLES  02  ABSORBED /KGMOLE 

1  i 

KB  = 

0  0490 

KCMOLES  N2  ABSORBED /KGMOLE 

1  1 

It  - 

2.6500 

1 

0 1  FBI  SION  COEFFICIENT  3 

200.00 

1/SEC 

.’  ID  -  RUCTION  3 

0.370 

t  J 

COB  RESPONDING  PARAMETER  #  TO  CHANCE  PARAMETER 

1  1 

■  "O'  TO  SIGNIFY  PARAMETERS  ARE  CORRECT— RUN  SIMULATION 

'  1 

1  "‘59"  -0  INSERT  A  STEP 

CHANGE  IN 

BREATHING  FLOWRATE 

1 

'  "J“"  0  EXIT  PROCRAM 

10 
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THE  FOLLOWING  PLOTS  ARE  AVAILABLE  TO  VOU  ON  THE  TERMINAL 


1  MOLE  FRACTION  OF  OXVGEN  VS  TIME  (DATA  F I LE=MOLEFR AC . DAT > 

2  INLET  MASS  FLOWRATE  VS.  TIME  (DATA  F I LE=FL0WR ATE . DAT > 

3  OUTLET  MASS  FLOWRATE  VS.  TIME  (DATA  FILE=FLOWRATE  DAT) 

4  A  DVNAMIC  SIMULATION  FOR  ONE  CVCLE  (OXVGEN  VS  TIME, 

DATA  FILE=PROFILE  DAT) 


PRESS  THE  CORRESPONDING  #  TO  HAVE  THE  OUTPUT  SHOWN  AS  THE 
SIMULATION  TRANSPIRES 


1 


1  ENTER  "0"  TO  RUN  SIMULATION 

2  ENTER  "19"  TO  CHANGE  THE  PARAMETERS 

3  ENTER  "99"  TO  EXIT  PROGRAM 


ENTER  THE  CORRESPONDING  *  AND  PRESS  "RETURN" 


0 


no  VOU  WISH  TO  DOUBLE  THE  NUMBER  OF  SPACE  LUMPS  IN 
! HE  SIMULATION  FOR  BETTER  ACCURACV  ?  (  CURRENT  NUMBER  USED 

in  50  )  (V/N) 
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SAMPLE  OUTPUTS 


«  OXYGEN  MOLE  FRACTION  « 

PSUP  *  16.69  PSIA 
POUT  »  1.69  PSIA 
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SAMPLE  OUTPUTS 


*  INLtT  HAS5  * 

PSUP  =  40.00  PSIA 
POUT  =  15.00  PSIA 


TIME  (SEC.) 


PHASE  PRESS  "C"  TO  CONTINUE  : 
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*  OUTLET  HASS  FLOWRATE  < KG/SEC)  « 


PSUP  *  40. 00  PSIA 
POUT  *  15.00  PSIA 
WBRL  *  10.00  STD  LIT/HIN 


TIME  (SEC.) 


PLEASE  PRESS  "C"  TO  CONTINUE  : 


OBOGS  SIMULATION 
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